Phonon bottleneck in graphene-based Josephson junctions at millikelvin temperatures 
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We examine the nature of the transitions between the normal and the superconducting branches of 
superconductor-graphene-superconductor Josephson junctions. We attribute the hysteresis between 
the switching (superconducting to normal) and retrapping (normal to superconducting) transitions 
to electron overheating. In particular, we demonstrate that the retrapping current corresponds 
to the critical current at an elevated temperature, where the heating is caused by the retrapping 
current itself. The superconducting gap in the leads suppresses the hot electron outflow, allowing us 
to further study electron thermalization by phonons at low temperatures [T < IK). The relationship 



between the applied power and the electron temperature was found to be P oc T 
is consistent with cooling due to electron-phonon interactions. 

PACS numbers: 74.45.+C, 72.80.Vp, 63.22.Rc, 65.80.Ck 



which we argue 



The electron-phonon interaction and the thermal prop- 
erties of graphene have attracted a lot of attention (for an 
extensive review, see e.g. Ref. [l|.) Most studies focused 
on relatively high temperatures; however, this year sev- 
eral experiments were performed at temperatures below 
10 Kelvin^ri. Measuring the intrinsic thermal properties 
of graphene becomes more challenging in this regime. In 
particular, making normal metal contacts to graphene 
would provide the dominant thermalization path, effec- 
tively shunting the electron cooling by phonons^. In or- 
der isolate the phonon contribution to thermal transport 
in graphene below IK, we have contacted the graphene 
crystal with electrodes made from lead (Pb), which be- 
comes superconducting below T ^ 7K. Superconducting 
leads exponentially suppress the thermal transport of hot 
electrons outside of the sample, allowing one to study 
thermalization by phonons. 

We study several supcrconductor-graphene- 
superconductor (SGS) junctions, and focus on the 
hysteresis between the switching current Is (from the 
superconducting to the normal state) and the retrapping 
current Ie. (from the normal to the superconducting 
state), which is observed despite the fact that our 
junctions arc overdamped^. Wc demonstrate that the 
hysteresis is due to electron overheating^iiS. Multiple 
contacts are made to the same graphene crystal, and re- 
sistive heating is applied at one location; simultaneously, 
we measure the electron temperature by monitoring the 
change in the critical current through a pair of contacts 
at a different location on the crystal. By changing the 
pairs of contacts where heating is applied and where 
the critical current is measured, we conclude that the 
graphene crystal is well thermalized, i.e. its electrons 
maintain a uniform temperature. We evaluate the elec- 
tron temperature as a function of applied power, which 
allows us to estimate the electron-phonon scattering rate 
in graphene at sub-Kelvin temperatures. 

Graphene was deposited on the Si/Si02 substrate us- 
ing the standard mechanical exfoliation recipe^^ and ver- 
ified to be a single layer by Raman spectroscopyi^. Six 
parallel superconducting strips (500 nm wide) were de- 



posited on top of the graphene crystal by thermally 
evaporating a 4 nm contact layer of palladium (Pd) fol- 
lowed by 120 nm of lead (Pb)i^. Thereby, a total of 
5 superconductor-graphene-superconductor (SGS) junc- 
tions were created with lengths of 0.3, 0.7, 1, 1.5 and 
2 /.im, labeled 1 through 5, respectively (Fig. la, b). 
The widths were roughly similar for all the junctions at 
around ^ 5fini. In order to measure the relatively small 
superconducting currents, the sample was encased in a 
copper box thermalized at the millikelvin temperature. 
The wires connecting the sample to the measurement 
setup were filtered by cold RC filters and thermocoax 
cables. 

At the base temperature of 35 mK the shorter junc- 
tions support a gate-dependent supercurrent (Figure Ic). 
The switching currents from the superconducting to nor- 
mal state Is at the fully open gate voltage of Vgate = 40 
V were ^ 40 nA, ^ 20 nA, and ^ 5 nA for junctions 
1 — 3 respectively. The two longest junctions (junctions 4 
and 5) never developed a supercurrent. These values are 
suppressed compared to the naive estimates of the criti- 
cal current Ic ^ A./Ri^e (Rn is the normal resistance of 
the junction) most likely due to disorder, as discussed in 
our earlier publicationi^. 

Nanoscale hybrid Josephson junctions tend to be un- 
derdamped due to the large capacitance shunting the 
junction's bonding pads to the back gate. (The damping 
factor Q of the Josephson junction is related to the criti- 
cal current Ic, junction capacitance C and the shunting 
resistance R as Q = R^y2j:IcG^.) We made the leads 
connecting the pads to the junctions somewhat resistive 
(a few hundred fl) to partially isolate the junctions from 
the pad capacitance and to introduce dissipation. As 
a result, we estimate that our junctions are in the over- 
damped regimei^. Indeed, successive current sweeps pro- 
duce a very narrow distribution of the switching and re- 
trapping currents (see inset to Figure 2b; also note that 
the variations of the switching current seen in Figure Ic 
are reproducible on successive sweeps of Vgate-) There- 
fore, the observed hysteresis between the switching and 
retrapping current cannot be due to an undcrdampcd dy- 
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FIG. 1. a) Scanning Electron Micrograph of the sample. A 
layer of graphene (light-colored region at the top of the im- 
age) is contacted by superconducting electrodes (dark verti- 
cal stripes) made from Pd/Pb bilayer (partially oxidized fol- 
lowing the measurements). Six contacts form five Josephson 
junction of different lengths (distances between the contacts). 
The superconducting electrodes extend just past one of the 
sides of the graphene crystal, where they are contacted by the 
normal metal electrodes. On their other end, the supercon- 
ducting electrodes do not completely cross the crystal, allow- 
ing a thermalization path to connect the different regions of 
graphene. b) Schematic of the measurement setup for Figure 
2. Differential resistance of one junction is measured, while a 
heating current Ih is applied to a different junction, c) Dif- 
ferential resistance map of junction 1 measured at the base 
temperature of 35 mK as a function of the back gate voltage. 
The supercurrent is visible as the central dark region. No- 
tice the hysteresis between the switching (positive) and the 
retrapping (negative) currents. 



namics, and we should seek an alternative explanation. 

In order to investigate the nature of the hysteresis, the 
graphene flake was heated locally by sending a current 
Ih through one of the five junctions, while measuring 
the critical current of another junction^. For example, 
junction 1 would be measured while the heating current 
Ih is sent through junction 5 (Figure 2a). At zero heating 
current, the values of the switching and retrapping cur- 
rents Is and Ir are the largest, and both decrease as the 
heating power Ph is increased. The difference between 
Is and Iji also decreases and the hysteresis disappears 
after some intermediate value of Ph ■ 

While the switching current Is is reached when the 
measured junction is in the superconducting state, the re- 
trapping current Ir is realized when the measured junc- 
tion is in the normal state. Therefore, at the retrap- 
ping transition, additional power Pr = RnIr is ap- 
plied to graphene. Taking into account this additional 
power, we plot Is{Ph) and Ir{Ph + Pr) in Figure 2b. 
(Here, Ph = RhIh is the power applied at the heater 
junctioniS,.) Ig and Ir clearly fall on the same curve. 
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FIG. 2. a) The switching and the retrapping currents of the 
shortest junction (#1) versus the heating current applied to 
the longest junction (#5). The difference between the switch- 
ing and the retrapping current disappears as more heating 
power is applied, b) Switching and retrapping currents ver- 
sus the total heating power at different gate voltages. Since 
the retrapping current is necessarily measured while the junc- 
tion is in the normal state, the power dissipated in that junc- 
tion contributes to the total heating power. With this ad- 
ditional heating taken into account, the switching and the 
retrapping currents fall on top of each other. We conclude 
that the hysteresis between the switching and retrapping cur- 
rents is caused by the self heating of the junction^ii^. Inset: 
successive current sweeps (negative to positive) of junction 2 
at Vgate ~ 40 V. The fluctuations in values of 1$ and Ir are 
negligible, indicating that the hysteresis is of a different origin 
compared to underdamped junctions. 



demonstrating that at the same total dissipated power. 
Is and Ir coincide. 

The curves of Is and Ir versus P map onto each other 
regardless of the junction that was measured (junctions 1 
or 2) and for different gate voltages applied to graphene. 
Therefore, we can conclude that the retrapping current 
is suppressed compared to the switching current due to 
the heat dissipated by the retrapping current itself; the 
hysteresis between the switching and retrapping currents 
is thus due to the self-heating of the measured junc- 
tion. Indeed, it has been demonstrated that self-heating 
is the origin of the hysteretic behavior in overdamped 
SNS junctions^ii^. This conclusion allows us to charac- 
terize the electron cooling in graphene in the remainder 
of the paper. 

Instead of heating the graphene crystal locally by the 
current Ih, we can heat the entire sample. We measure 
the switching current Is in the two shortest junctions 
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FIG. 3. Dissipated power vs. electron temperature in 
graphene. The latter was extracted by matching the switch- 
ing current at a given heating current to the switching current 
measured by heating the whole sample. The graph presents 
three measurement configurations: a) measuring the switch- 
ing current of junction 1 while heating junction 3 or junction 
5, and b) measuring junction 2 while heating junction 1. The 
data for all three configurations fit well to a P oc T"^ power law 
dependence with less than 5% difference in the proportional- 
ity factor, meaning that the electrons in the entire region are 
thermalized. 



versus temperature T at several difFercnt gate voltages. 
By comparing the Is vs. T data with the I3 vs. P data, 
we are able to extract the temperature up to which the 
electrons in the measured junction are heated for a given 
P. This way, a plot of P vs. T has been extracted for 
measured junction 1 while heating either junction 3 or 5 
(Figure 3a). Similarly, in Fig 3b, we measured junction 2 
while heating junction 1. (In the latter case, large enough 
current was passed through junction 1 to overcome its 
own supercurrent. Once that happens, graphene heats 
to about 200 mK, which explains the gap in the data at 
lower temperatures.) 

The three curves in Figure 3 are remarkably close, al- 
though the three setups have different distances and con- 
figurations of the heater and the measured junctions, im- 
plying that the entire region of graphene crystal in the 
vicinity of the junctions is thermalized to a uniform elec- 
tron temperature. Empirically, the curves can be fitted 
by the same power law P = 6 • lO^^^ v/ rj.^^^ 
pendence rules out cooling by the hot electron diffusion 
into the superconducting leads. Indeed, the dependence 
would be exponentially suppressed in temperature due 
to the gap in the quasiparticle density of states of the 
superconductor. Even at T =1 K, the dissipated power 
(6 ■ 10~^^ W) is several orders of magnitude larger than 
that estimated from the Wiedemann-Franz law, properly 
adjusted for the suppressed density of states in the leads 
(^ 10^^^ W). This situation is different from the regime 
realized in Ref. [H, where a relatively small superconduct- 
ing gap in Al leads allows them to provide a dominant 
thermalization path. We conclude that in our sample 
electrons in the graphene crystal are thermally decoupled 



from the leads and must be cooled by phonons. 

Furthermore, the observed P{T) cannot be limited by 
the Kapitza resistance between graphene and the sub- 
strate, which is estimated to be a few orders of magni- 
tude smaller. (As an upper limit for the resistance, we 
take the values measured at 40 K in Ref. [l^ and extend 
the cooling power according to cx or oc T*.) Hence, 
graphene lattice must be well thermalized with the sub- 
strate. We conclude that the bottleneck for cooling in 
our system is the weak thermal coupling between elec- 
trons and phonons in graphene. 

Following the derivation of Ref. [13, we estimate 
the electron-phonon cooling power in graphene as 

h^p^, i O^bTeiY . Here, A is the area of graphene, D 
is the deformation potential constant, n is the electron 
density, p,n is graphene mass density, s is the speed of 
sound, and Tg/ is the electron temperature. We assume 
that at low temperatures the wavelength of the emitted 
phonon hs/ksTei, which enters the scattering matrix el- 
ement, should be replaced by a distance I of the order 
of a hundred nm {e.g. the distance between the leads, 
or the electron mean free path). We also took into ac- 
count that Ref. [l3 calculated the transport time while 
we are interested in the scattering time, because each 
phonon emission event typically results in the electron 
energy loss of the order of ksTei. 

Plugging in the numbers from Refs. [13 and[l^ we es- 
timate the cooling power at P = 10~^^^ T^^, within an 
order of magnitude from the measured value. [This is a 
rather crude estimate, in part because we do not know 
how large is the thermalized area of graphene crystal, 
which extends beyond the immediate vicinity of the con- 
tacts, see Figure la. This area determines A, taken to 
be (10/xm)^ in this estimate.] 

Recent measurements of Refs. and report cool- 
ing power for hot electron P cx T^i, indicating that the 
emitted phonon wavelength hs/ksTei is not cut off, as 
we assumed above. The difference may stem both from 
the lower temperatures in our measurement and a more 
restricted sample geometry resulting in shorter I in our 
case. The cooling power of Ref. 0, extrapolated down to 
1 K would yield 0.07 W/m^, which in fact is close to our 
result for A = (lO/im)^. 

Finally, P cx T^i dependence due to supercoUision 
coolingi^ has been just observed in Refs. [l^ and [21I . 
This regime is realized if the wavevector of the emitted 
phonon, fc^Tgi/s, exceeds the electron Fermi wavevector 
kp; the emission process is enabled by the disorder. We 
work in the opposite regime, kp ^ ksTei/ s, so that the 
theory of Ref. [l9| cannot be directly applied. 

In conclusion, we have shown that the difference be- 
tween the switching and the retrapping currents in our 
graphene Josephson junctions is caused by electron over- 
heating in the normal state. The superconducting con- 
tacts thermally isolate the graphene crystal from the 
leads, allowing us to measure the electron temperature 
rise for a given dissipation power, and hence the electron- 
phonon energy transfer rate. The observed power law 
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dependence P oc is consistent with theory of clectron- 
phonon interactions in graphene and with other measure- 



ments. Due to small electron heat capacitance and their 
decoupling from phonons, this type of sample may be 
useful for detector applications^^. 
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